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ABSTRACT 

We present an analysis of the proximity effect in a sample of ten 2 A resolution QSO 
spectra of the Lya forest at (z) = 2.9. Rather than investigating variations in the 
number density of individual absorption lines we employ a novel technique that is based 
on the statistics of the transmitted flux itself. We confirm the existence of the proximity 
effect at the > 99 per cent confidence level. We derive a value for the mean intensity 
of the extragalactic background radiation at the Lyman limit of J = S.Sl^'g x 10~^^ 
ergs s""'^ cm~^ Hz~^ sr~^. This value assumes that QSO redshifts measured from high 
ionization lines differ from the true systemic redshifts by Av ^ 800 km s^^. We find 
evidence at a level of 2.6cr that the significance of the proximity effect is correlated with 
QSO Lyman limit luminosity. Allowing for known QSO variability the significance of 
the correlation reduces to 2.1cr. 

The QSOs form a close group on the sky and the sample is thus well suited for an 
investigation of the foreground proximity effect, where the Lya forest of a background 
QSO is influenced by the UV radiation from a nearby foreground QSO. From the 
complete sample we find no evidence for the existence of this effect, implying either 
that J > 20 X 10~^^ ergs s~^ cm~^ Hz~^ sr~^ or that QSOs emit at least a factor of 1.4 
less ionizing radiation in the plane of the sky than along the line of sight to Earth. We 
do however find one counter-example. Our sample includes the fortunate constellation 
of a foreground QSO surrounded by four nearby background QSOs. These four spectra 
all show underdense absorption within ±3000 km s~^ of the redshift of the foreground 
QSO. 
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1 INTRODUCTION 

The study of many physical processes at high redshift re- 
quires knowledge of the intensity of the UV background ra- 
diation, J. For example, it is thought that the Lya forest in 
QSO absorption spectra is caused by highly photo-ionized 
gas and thus an estimate of the total mas s content of the in- 
tergalactic medium (I GM) depends on J (Ranch et al. 1997 
Weinberg et al. 1997). It is also one of the paramet ers that 
define the environment in which galaxies form (e.g. 



Susa & 



Umemura 200C), and its value and evolution provide impor- 



tant constraints on t he objects believed to be the origin of 



the background (e.g. Bechtold et al. 1987 Haardt & Madau 



199e; Fardal, Giroux, fc Shull 199^). At high redshift, the 



background is often measured from the proximity effect, i.e. 
the observed underdensity of Lya forest absorption lines in 

the vicinity of backgroun d QSOs. 

Nearly 20 years ago Carswell et al. (1982| ) first noted 



that the mean density of Lya absorption lines seemed to 
increase with redshift when comparing the spectra of several 
differ ent QSOs and yet dec rease along individual lines of 
sight. Murdoch et al. (1986) confirmed this 'inverse effect', 
established that it was confined to the vicinity of the QSO 
('proximity effect') and offered two possible explanations: 
i) the absorbers near a QSO may be too small to fully cover 
the continuum emitting region or ii) absorbers in the vicinity 
of a QSO may be more highly ionized than elsewhere due to 
the QSO's UV radiation. 



In a seminal paper Bajtlik, Duncan, & Ostriker (198J 



E-mail: jol@phys.unsw.edu.au 



(hereafter BDO) developed a quantitative ionization model 
and, for the first time, measured J from the observed un- 
derdense absorption near 19 QSOs and their observed lu- 
minosities. In the most comprehensiv e intermediate resolu- 
tion study to date Scott et al. (2000[ ) analysed a sample of 
74 spectra. Like Bechtold (1994 ), they divided their sample 
into low and high-luminosity subsamples and found that the 
relative deficit of absorption lines within 1.5 h^^ Mpc of the 
background QSOs was more significant in the latter. 
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Figure 1. Summary of previous proximity effect measurements 
of tlie mean background Lvman limit intensity. J. in units of 
10~^^ e l'gs s~^ cm~^ Hz~^ sr~^. Solid and open symbols rcp- 
resent measurements from intermediate and high resolution data 
respectively. Estimates from the background (foreground) prox- 
imity effect are shown as circles (triangles), those that include 
both are shown as squares. Error bars are la. The dashed line 
marks the value of J at 2 = 2.5—3 computed by Haardt & Madau 
(1996) for a background dominated by the observed QSO pop- 
ulation and for a go = 0.1 cosmology. The top panel shows the 
redshift intervals covered by the various studies. 



CoDke, Espey, & Carswell (1997) studied the proximity 



effect at high spectral resolution. In their thorough analysis 
of 11 spectra they gave a detailed account of the various 
statistical and systematic errors and biases involved in the 
measurement of J from the proximity effect. Like previous 
authors they found no evidence to suggest that the back- 
ground intensity evolved with redshift. 

We summarize these and other measurements of J in 
Fig. 1^. It is worth noting that not all the points in this plot 
are independent of one another. There are large overlaps in 
the data used. For example, Q0014-)-813 is included in eight 
of these studies. In addition, all except one of these mea- 
surements are based on 'line cou nting', i.e. the statistics of 



individual absorption lines. Only Lu et al. (1996 ) considered 
the integrated absorption in 100 A bins. [Zuo (1992 ) obtained 
a rough estimate of J from correlations (where W i s 

the rest equivalent width) and M0ller & Kjargaard (1992) 



used a basic flux statistics approach to investigate the fore- 
ground proximity effect. Thus no alternative methods seem 
to have been explored in great detail. 

Most authors have found good agreement between their 
data and the ionization model of BDO and alternative ex- 



planations for the proximity effect have not received much 
observational support. The sizes of Lya ab sorbers inferred 



from observations of close QSO pairs (e.g. Dinshaw et al 



199S and references therein) seem to rule out the possibil 
ity that the absorbers are too small to completely cover th e 
background QSO. In addition, Lu, Wolfe, & Turnshek (1991) 



found no difference between the W distribution of lines near 
QSOs and that of lines far from QSOs. They also eliminated 
a broken power law for the redshift distribution of lines as 
a possible cause for the proximity effect. 

Thus increased ionization due to the extra UV flux from 
the QSO seems to remain as the only credible explanation 
for the proximity effect. However, it implies two observable 
effects: 

1. The proximity effect should correlate with QSO lu- 
minosity. More luminous QSOs should deplete larger regions 
more thoroughly than less luminous ones. BDO claimed that 
this effect was present in their data and that it was con- 
sistent with the exp ectations from the ir ionization model. 
Bechtold (1994|) and IScott et al. (2000|) also found a weak 
correlation. Lu, Wolfe, & Turnshek (1991) on the other hand 
found no evidence for a correlation at all but nevertheless 
concluded on the basis of simulations that this was consis- 



tent with the ionization model, srianand & Khare (199(:) 
could not identify a correlation either. 

2. In addition to the 'classical' background proximity ef- 
fect there should be a foreground proximity effect where the 
absorbing gas along the line of sight to a background QSO is 
depleted by the UV radiation of a close-by foregroun d QSO. 
Study ing a triplet of QSOs separated by 2 to 3 arcmin [Srotts 



found no evidence for the existence of the foreground 



proximity effect. M0ller & Kjsrgaard (1992) added another 



spectrum to this triplet and conflrmed the negative result. 
Dobrzycki & Bechtold (1991) observed a ~ 10 Mpc void in 
the spectrum of Q0302-003 with a foreground QSO sepa- 
rated from its line of sight by 17 arcmin. However, the fore- 
ground QSO was displaced from the void by ~ 3600 km s~^, 
implying either that the QSO radiates anisotropically or 
that it turned on on a time-scale co mparable to the light 



travel time from the QSO to the void. Fernandez- Soto et al 
(1995| ) studied three QSO pairs separated by 3.8 to 12.6 ar 



cmin but they were unable to rejec t the non-existe nce of the 
effect by more than ~ 1 a. Finally, srianand (1997) reported 
a ~ 7 Mpc void in the spectrum of Tol 1038-2712 with a 
foreground QSO at t he redshift of the void and separated 
by 17.9 arcmin. Like Dobrzycki & Bechtold he showed that 



it was unlikely that the void was a chance occurrence. Thus 
there currently exists only a single example where an under- 
density of absorption lines in the spectrum of a background 
QSO can be explained by the presence of a foreground QSO 
without making extra assumptions. 

The main goals of this paper are to introduce a new 
method to analyse QSO spectra for the proximity effect and 
to address the two problems described above. The data used 
in this investigation are described in Section ^. They consist 
of the spectra of a close group of 10 QSOs which have not 
been included in any previous studies. Thus the analysis 
presented here is independent from others in the sense that 
it uses both a different method as well as different data. In 
Section ^ we measure J from the classical proximity effect 
and demonstrate its correlation with QSO luminosity. In 
Section ^ we turn to the foreground proximity effect. Finally, 
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Table 1. QSOs. 



RA offset (arcmin) 



Figure 2. Distribution of QSOs in the sky. Tlio field is centered 
ona = 00h42™10° and (5 = -26°40' CB1950). Solid stars marltthe 



positions of the QSOs whose spectra are analysed in this paper. 



Empty stars mark the positions of additional foreground QSOs 
in the field. Emission redshifts are indicated. 



we consider a range of uncertainties in Section ^ and discuss 
our results in Section ^. 

Unless explicitly stated otherwise we use 50 = 0.5, A = 
and Ho = 100 h km s~^ Mpc"^ throughout this paper. 



2 THE DATA 



Williger et al. (1996) performed a large Civ survey in the 
South Galactic Pole region. Here, we use that subset of 
the data which has reasonable coverage of the Lya forest 
(10 spectra). The instrumental resolution is ~ 2A and the 
signal-to-noise ratio per pixel reaches up to 40 per 1 A pixel. 
The obse rvations and the redu ction process are described in 
detail by [Williger et al. (1996|). 



A search of the literature revealed seven additional 
QSOs in this field and in the appropriate redshift range. 
These will be considered as potential foreground ionizing 
sources. The angular separations range from 6.1 to 95.8 ar- 



cmin atid the emission redshuts range trom to ^.44. 
The di^Li'ibutlou of aU of Lheae QSOy in tho ijky in Aiawii lu 



Fig-i 



Since absolute spectrophotometry was not available for 
any of these QSOs we had to estimate continuum ilux den- 



Object 






a 


L^(912)'' 


Refs 




Onn41 -26(17 


2.505 


17.23 


—0.60 


15.01 


1 2 






3.053 


18.35 


—0.53 


12.17 


1 2 




00041 -2fi^S 


2.457 


18.70 


—0.88 


2.34 


1 2 




00041—9707 


2.786 


18.03 




10.14 


1 2 




00042-2627 


3.289 


18.55 




13.29 


1 2 




O0042— 9fi'^Q 


2.98 


20.05 


—0.09 


4.59 


3 




Q0042-2656 


3.33 


19.55 


-0.80 


4.11 


3 




Q0042-2657 


2.898 


18.78 




5.88 


1, 2 




Q0042-2714 


2.36 


19.88 




1.09 


4, 2, 


5 


Q0043-2633 


3.44 


19.61 


-0.77 


5.01 


3 




Additional fore 


ground 


QSOs: 










Q0040-2606 


2.47 


19.48 




1.80 


4, 2 




Q0041-2622 


2.18 


19.28 




1.69 


4, 2, 


5 


Q0042-2642 


2.81 


20.60 


-1.40 


0.28 


3 




Q0043-2555 


3.31 


20.53 


-0.89 


1.43 


3 




Q0043-2606 


3.11 


20.37 


-0.96 


1.07 


3 




Q0044-2628 


2.47 


19.28 




2.18 


4, 2, 


5 


Q0044-2721 


3.16 


20.17 


-0.45 


3.05 


3 





"The typical error on ms is 0.15. 
''In units of 10'^'' ergs s~^ Hz 
References: (1) R edshifts and Bj magnitudes from 



. _ ^ Hewitt, Foltz,, 

Chaffee (1995| ). Where given, continuum slopes were measured 



from low-resolution spectra kindly provided by Paul Francis. 
(2) We converted from Bj to Johnson B using the colour equa- 
tion oflBlair fc Gilmore (198S) and as suming {B — V) = 0.3. 



(3) Warren, Hewitt, & Osmer Q991 ). Conversion to i? magni 



tud cs using the colou r equation of 
(4) [Prinkwater (1987|). 



Warren et al. (1991) 



(5) Redshifts from Williger et al. (1996) 



-0. 6 this equation gives a flux 1.6 times highe r than 



Tytlei's (1987) empirical formula. B[owever, note that Tytler 



used the AT-corrections of Evans fc B[art (1977 ) whereas we 
use the if-corrections given by [Cristiani fc Vio (199C ) be- 
cause they extend beyond z = 2.5. For some QSOs con- 
tinuum slopes were not available. In these cases we used 
a — —0.6 which is similar to the value of —0.5 given by 



Francis (1993| ). 

We correct the above flux value for Galactic extinction 
by applying a correction factor lo" *'*'^' where 



A{X) = Rv E{B - V) 



A(\) 

A{vy 



(2) 



We use Rv = A{V)/E{B - V) = 3.1 w hich is the average 



value for the diffuse interstellar medium ( Clayton & CardeUi 



1985). The variation of the extinction wit h wavelength rela - 



tive to that a t V, AiX)/A(V), is given by P'Donnell (1994| ) 
(optical) and [Cardelh, Clayton, fc Mathis (1989D ( UV). We 
take E{B — V) for each QS O from the dust map of schlegel 



sities fruiii ubstifved i3-baiid maguUudys. AsHUiiiiiig a power 
law continuum oc i^" the observed flux density fi, at ob- 
served wavelength A is given by 



Finkbeiner, fc Davis (1998|) . This procedure is equivalent to 



/4A) = 



A, 



A(H-zq)-i 



(1 + zq)10-°-*('"^-'=^V.x(0)(1) 



where Ax, mx, kx and f,yx{0) are the central wavelength, 
observed m agnitude, if-correction and 0-magnitude flux 
(Allen 1991) of the X-band respectively. For X = V and 



first correcting m_B for Galactic extinction and then apply- 
ing a correction to a due to the wavelength dependence of 
the extinction. 

Gaskell (1982) first pointed out that QSO redshifts mea- 
sured from high ionization emission lines like Lya or C iv are 
systematically lower than redshifts measured from lower ion- 
ization lines like Mg 11 or the Balmer series which are thought 
to indicate the systemic redshifts. When estimating J from 
the classical proximity effect using QSO redshifts derived 
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from high ionization lines, the resuh will be too high because 
the lower QSO redshift implies a higher QSO flux at a given 
cloud and th erefore (for th e same observed effect) a higher 
background. Espey (1993) showed that a velocity shift of 



1500 km s lowered the value of J from log J — —20.75 



Lu, Wolfe, fc Turnshek's (1991) data 



to -21.30 in 

The redshifts of the QSOs considered in this paper were 
all determined from high ionization lines. However, since 
there is considerable disagreement in the literature over the 
values of the line shifts and possible correlations with QSO 
luminosity and/or emission line properties (e.g. Tytler & 
Fan 19^ and references therein), it is difflcult to reliably 
correct for this effect. We shall therefore resort to determin- 



ing J as a function of line shift in Section 4.2. 

In Table ^ we list all the QSOs considered as well as 
their redshifts, _B-band magnitudes, continuum slopes and 
Lyman limit luminosities calculated from equations (|l|) and 

&■ 



3 ANALYSIS 
3.1 Method 



Essen tially, we use the method of Liske, Webb, & Carswell 



(1998) (hereafter LWC) to look for regions of underdense ab- 



sorption (i.e. 'voids') in Lya forest spectra near the positions 
of suspected sources of ionizing radiation. These sources may 
be the background QSOs themselves or foreground sources 
near the line of sight to the background QSOs. The tech- 
nique does not rely on investigating variations of the num- 
ber density of individual absorption lines but rather uses the 
statistics of the transmitted flux directly. Thus it elegantly 
sidesteps all problems related to the incompleteness of lines 
due to limited spectral resolutio n and signal-to-noise, W - 



limite d versus A'^-limited samp les ( Chernomordik fc Ozernoy 



199; ; Srianand fc Khare 1996) and Malmquist bias (Cooke 



Espey, fc CarsweU 1997). M0ller fc Kjaargaard (1992) have 



previously employed a somewhat more rudimentary version 
of our method to investigate the foreground proximity ef- 
fect. Flux statistic s have also been used in othe r contexts b y 
Liske ct al. (2000| ). |Nusser fc Hachnch (1999|). |Hui (1999|) , 



Croft ct al. (19991), [Zuo fc Bond (199 j|). |Zuo fc Lu (1993 ) 



Press, Rybicki, fc Schneider (1993|), [Press fcRybicki (1993D 



Webb et al. (1992] ) and [Jenkins fc Ostriker (1991] ). 

In order to identify large-scale regions of underdense ab- 
sorption we convolve a spectrum with a smoothing function 
in order to filter out the high-frequency 'noise' of individ- 
ual absorption lines. A large-scale feature will be enhanced 
most by this convolution if the scale of the smoothing func- 
tion, (Ts, is similar to the scale of the feature. Since we do 
not know the scale of any possible features a priori we sim- 
ply smooth the spectrum on all possible scales, from 1 pixel 
where the data remain unchanged to the largest possible 
scale where the data are compressed into a single number 
akin to 1 - Da ( |Oke fc Korycansky 198^ ). When plotted 
in the wavelength-smoothing scale plane this procedure re- 
sults in the 'transmission triangle' of the spectrum which 
we denote by G(A, as) for the case of a Gaussian smoothing 
function and by T(A, cts) for the case of a top-hat smoothing 
function (Section 3.4). 

We then compare the observed transmission to the ex- 



pected mean transmission calculated on the basis of the sim- 
ple null-hypothesis that any Lya forest spectrum can be de- 
scribed as a collection of absorption lines whose parameters 
are uncorrelated. We also calculate the expected variance of 
the transmission in order to assess the statistical significance 
of any fiuctuations of the observed transmission around the 
expected mean. 



3.2 Basic absorption model 

LCW considered a (normalized) Lya forest spectrum as a 
stochastic process where each point in the spectrum is a ran- 
dom variable, e~^, drawn from the transmission probability 
density function. They calculated the mean and variance of 
this function on the basis of the above null-hypothesis. In- 
cluding instrumental effects in their calculation they then 
used these results to derive the mean and variance of G, 



{G)(A,a, 



and 



crG(A,crs) 



where 



cj^-r = exp 



exp 



-B 



^n(A) 

20r as/pf 



+ 



^ 5^ (A) + ^ 



-Hi) 



exp 



-2B 



7 + 1 



(3) 



(4) 



.(5) 



(j„ denotes the noise of a spectrum, ctlsf the width of the 
instrumental line spread function, ps the pixel size in A 
and Xa — 1216 A. q is the intrinsic width of the auto- 
covariance function of a 'perfect' spectrum (i.e. before it 
passes through the instrument) whic h we determine from 
simulations (LWC; [Liske et al. 2000| ). 



B contains the nor- 



malization and is given by 



N'^f{h) W{N,b) ANdb, 



(6) 



where W[N, b) is the rest equivalent width of an absorption 
line of column density A'^ and Doppler parameter b. B will 
be measured directly from the data (Section 3.4). A, 7,/?, 



and }{b) are all part of the observationally determined dis- 
tribution of absorption line parameters 



AzANAb 



A{l + zyN-''f{b). 



(7) 



We take f rom recent high resolution studie s as 1.5 ( Hu 



et al. 1995[; Lu et al. 1996, Kim et al. 1997; Kirkman fc 



Tytler 1997| ). For the present data set [WiUiger et al. (200C| ) 
determined 7 = 2.0 which we will adopt. We di scuss the 



effects of uncertainties in /3 and 7 in Section 4.2.1 



3.3 Incorporating local ionizing sources 

In order to measure J it is necessary to extend our model 
of the transmitted fiux to incorporate local fiuctuations of 
the ionizing radiation caused by discrete sources. We adopt 
the simple ionization model of BDO which basically consists 
of the assumption that the absorbing gas is highly photo- 
ionized so that an absorber's column density is inversely 



The proximity effect in a close group of QSOs 5 



proportional to the incident ionizing flux. Thus in regions of 
enhanced ionizing radiation we must modify equation (Q): 



d^AA 



oc {l + z)'<N-'^{l + Lj{z)]^- 



(8) 



(9) 



d^diV 
where 

'^^ 4ttJ{z) 

J^i^(All) is the flux from the ionizing source (IS) received by 
the absorber at wavelength All = 912 A in the restframe of 
the absorber. The IS may be the background QSO itself or it 
may be a different, foreground QSO. The validity of equation 
(^) is subject to the limitation that the spectral shape of the 
background J below the Lyman limit is similar to that of 
the IS. This will be the case if the QSOs are the dominant 
cont ributors to t he background, if the IGM is optically thin 
(see Espey 1993 for a discussion of the optically thick case) 
and if the emission from the I GM does not drastically alter 
the shape of the background (Haardt & Madau 1996). 
Fi, may be calculated from 



F^All) = 



47r rl(za,z{^) 



(l + ^Is). 



(10) 



r-L{za,Zig) denotes the luminosity distance between the ab- 
sorber and the IS which, in general, is a function of the ab- 
sorber redshift, Za, and the redshift of the IS as seen by the 
absorber, z^g. Since we will consider the foreground prox- 
imity effect, where the absorber does not lie along the line 
of sight to the IS, zfs is in turn a function of Za,zis, and 
the angle a by which the abs orber and the IS are separated 
on the sky. See Liske (2000) on how to calculate z{g and 
rh{za, zis,a). Note the bandwidth correction factor (1 -l-ajg) 
which is usually ignored at this point. The intrinsic lumi- 
nosity of the IS, L^, is related to the observed flux at the 
observed wavelength A by 



V 1 + 2IS / 



4n rl(zis) 
1 + ZlS 



(11) 



where rL{zis) is the luminosity distance from Earth to the 
IS. Thus we have 



Uj{Za, zis, a) = 



/.(a 



ATTj{Za) 1 -I- ZIS 



r-bizis) 



rL{za,zis,a) 



.(12) 



The above ionization model is easily incorporated into 
our transmission model. Because of equation m the modi- 
fication M) implies 



B B{Za) = B [1+ Lu{Za, zis , a] 



11-13 



(13) 



It is unlikely that we will be able to constrain the red- 
shift evolution of the background with the present sample as 
previous studies of similar size but larger redshift coverage 
have been unable to do so (Cooke, Espey, & Carswell 1997 
Giallongo et al. 1996). However, the background is expected 



to peak smoothly in the redshift range covered here (Haardt 
& Mad m 199f: ) and the near constancy of J at z > 2 has 



been s upported bv the consistencv of simulations with the 
redshift ] evolution of the Lya forest ( Dave et al. 1999 ). For 



these reasons we assume J(z) — const. 



3.4 Further modifications to the absorption model 

The inclusion of the proximity effect in our model ren- 
ders one of the approximations of LWC in the derivation 
of equations (|^) and (^) invalid. There we used the fact 
that the mean (and the variance) of the transmitted flux, 
{e~^) = exp[— i3(A/Aa)^^^], is approximately linear in A 
over the scales of interest, so that, e.g.. 



(G)(A,a.) = 



2% o"s 



exp 



-B 



(e ^)(A') exp 
A N^+i' 



(A -Ay 
2a| 



(a,,) 



dA' 



(14) 



Because of the introduction of u>, approximations like the 
one above are no longer valid and we must now carry out 
all convolutions explicitly: 



{G)j(A,a.) = 



X exp 



27r (Jlsf (Ts 
(A' - A")^ 

2n-2 
^"LSF 



exp 



-B{z) ( f 



7+1 





r (A-A')^i 


exp 


2ai 



and 



o"Gj(A,o-s) = 



X exp 



\/2tv (Jc 
(A -AT 



2ai 



a^(A') 
dA', 



+ 



dA"dA' (15) 



-e'-.(A') 



(16) 



where as — a^l\f2. 

When investigating the classical proximity effect it will 
be helpful to use a top-hat smoothing function rather than 
a Gaussian because we are working at the 'edge' of the data. 
The equivalent of equations (|l^) and (^6|) are given by 



{r)j(A,as 



X exp 



\/27r ctlsf 2(Ts ja-o-b 
(A'- A")'' 



exp 



-B{z) 



7 + 1 



^"lsf 



dA"dA' 



(17) 



and 



o-|j(A,o-s) = - — 



A + tTa 



A — o-g 



a^(A') , ag-.(AO 
2crs/ps 



2as 



2ct3-2|A-A' 

exp 



2 [2a 



2 

LSF 



-g2(A' 



dA" 



dA'(18) 



As one approaches the background QSO in the classical 
proximity effect the flux from the QSO increases and B(z) 
decreases. Thus very close to the QSO the model predicts 
a mean transmission of almost 1 and a variance of almost 
0. This is clearly unphysical as absorption lines with Za ~ 
zq and even with Za <; zq are frequently observed. One 
of the reasons for this o bservation may be that absorbers 
have peculiar ve locities (Srianand & Khare 1996; Loeb & 



Eisenstein 1995 ). We accommodate peculiar velocities by 
convolving B{za) with a Gaussian of width 300 km s~^. 

We determined the normalization constant B directly 
from the data. First, we excluded all spectral regions with 
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c<j > 0.1 from background QSOs, assuming a fiducial value 
of log J = —21.0. For the present sample, the average size 
of the excluded regions translates to 5000 km s~^. For the 
remainder of each spectrum we then computed r((Js,inax) 
where (Ts^max is the largest possible smoothing scale. To these 
we then fitted our absorption model with B as a free param- 
eter. 



4 RESULTS ON CLASSICAL PROXIMITY 
EFFECT 

4.1 Significance 

With an absorption model and all its parameters in place 
we can proceed by transforming the transmission triangle 
of a given spectrum, T{\,aa), to a 'reduced transmission 
triangle' (RTT) by 



can already gain a useful estimate by simply maximising the 
expectation value of equation (EQ), 



RTj{\,as 



T-{T), 



(19) 



The reduced triangle has the mean redshift evolution of the 
absorption removed and shows the residual fiuctuations of 
the Lya transmission around its mean in terms of their sta- 
tistical significance. When neglecting the proximity effect in 
the calculation of (T) (i.e. J = oo), its presence in the data 
should be revealed by a region of -RToo > near the red 
edge of the triangle. 

We can consider the entire dataset at once in a com- 
pact manner by constructing a combined RTT: we first shift 
the spectra into the restframes of the QSOs, construct their 
transmission triangles and then average them where they 
overlap. Since different lines of sight are uncorrelated the 
variance of this composite is essentially just /n, where n 
is the number of spectra used. Thus at rest wavelength Al- 
and at restframe smoothing scale a^r we have 

1 ^ 

RTooiK, (Tsr) = — > {Ti[{l + ZQi)Xr, (1 -|- ZQi)asr] 

n ^ — ' 

1=1 

" (T) J = cx)[(l + ZQi)\r, (1 + ZQt)asr]} 

-1—1 

^ n 2 
-—yaT,J=oo[{i + ZQi)Xr,{l + ZQi)asr] , (20) 

where Ti and ZQi are the measured transmission and redshift 
of the ith QSO and we use J — oo in the calculation of (T) 
and (TT. 

The most significant 'void' in this combined RTT lies at 
(Ar,a"sr) = (1208.6 A, 3175 km s~^), two pixels away from 
the red edge of the triangle where one would expect a sig- 
nature from the proximity effect. It is significant at a level 
of 5.2(jt. Excluding each of the individual spectra from the 
composite in turn results in the significance of the feature 
varying from 4.2itt to 6.2(tt, so that the effect is not dom- 
inated by any single spectrum although there seems to be 
some variation in the strength of the effect among the in- 
dividual spectra. This will be investigated in more detail in 
Section 



4.3 



In this context it is helpful to ask what sort of signal 
one would expect from the data if the ionization model of 
Section 3.3 were correct. In the next section we will answer 



['^ Si = l ^TJ = oa] ^ 



(21) 



this question in detail with the help of simulations but one 



with respect to (Ar,o-sr). For example, for J21 — 1 (where 
J21 = Jx 10^^ ergs s~^ cm"'^ Hz~^ sr~^) we find a maximum 
of 5.0o"t. Note that the expected significance of the signal 
is not so much a function of the data quality but rather of 
the number of spectra included in the analysis, since the 
variance of the transmission is dominated by the 'noise' of 
individual absorption lines. 

An alternative explanation for the observed 'void' would 
be a systematic overestimation of the QSO redshifts. In this 
case we would expect to see an effect similar to the one ob- 
served because we would be including parts of the spectra in 
our analysis which correspond to regions physically behind 
the QSOs and thus would show much less absorption than 
expected. Note, however, that the redshifts of the QSOs were 
determined from high ionization lines and are thus expected 
to be too low, as discussed in Section ^, and not too high. 

We conclude that we have detected a proximity effect 
at a significance level of > 99 per cent. 



4.2 Measurement of J 

The seemingly most straightforward way to derive an es- 
timate of J would be to directly fit our absorption model 
to the observed composite transmission triangle with J as 
a free parameter. Fitting the transmission triangle rather 
than just the spectra has the advantage of ensuring that 
the model fits on all scales. In the previous section we saw 
that the data deviate most significantly from a no proximity 
effect model at a scale of ~ 3000 km s^^. On the smallest 
smoothing scale the strongest deviation is only 3o"t . Thus we 
can anticipate that by considering all scales we may derive 
tighter constraints on J. However, the pixels of a transmis- 
sion triangle are strongly correlated with one another and 
thus one would need to specify the entire covariance ma- 
trix in order to judge the quality of a fit. Instead we will 
use a much simpler yet effective approach which consists of 
considering only the most significant positive and negative 
deviations of the model from the data as a function of J. 

We implement this approach by searching for the most 
significant local extrema of an RTT in the region most likely 
to be affected by the proximity effect. A local maximum 
(minimum), LMAX (LMIN), is defined as any pixel in the 
RTT with RTj > {RTj < 0) and where all adjacent pix- 
els have smaller (larger) values. If there is more than one 
LMAX (LMIN) in a given wavelength bin (but at different 
smoothing scales) only the most significant one is consid- 
ered. 

We then define the search region as all those pixels in 
the RTT for which {RTao}{J) is larger than some threshold 
value ri. We begin by setting — I.Octt and search for the 
most significant LMAX and LMIN in this region. If none are 
found we increase the size of the search region by decreasing 
r] until the first LMAX (LMIN) has been found. Note that in 
order to calculate {RToo){J) we need to assume some value 
for J, which is precisely what we are trying to measure. 
However, starting with J21 = 1 the procedure converges 
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after only two iterations. Even without iterating, the above 
procedure ensures that the result does not depend sensitively 
on the exact value of J (within sensible limits) chosen to 
define the search region. 

In Fig. ^ (main panel) we plot both the most signifi- 
cant positive and negative deviations of the model from the 
data detected in this way as a function of the value of J 
that was used in the construction of the RTT. The void dis- 
cussed in the previous section is represented by the dot at 
J — oo (model with no proximity effect). For models with 
log J21 > 0.0 the RTT of the data shows significant un- 
derdense absorption (i.e. maxima at > 2(tt). On the other 
hand, for small values of J the model predicts too little ab- 
sorption on smaller scales and the RTT of the data shows 
significant overdense absorption (i.e. significant minima) for 
models with log J21 < 0.1. Thus there seems to be no value 
of J for which the model is entirely consistent with the data. 

However, recall that we are considering the maximally 
deviant points. For the correct model, the expectation value 
of the difference between a randomly chosen data point and 
the model is 0. However, given the additional information of 
the data point's rank, this is no longer true. For example, 
the probability distribution function of the maximum of a 
set of n uncorrelated, normally distributed numbers is given 
by P{x) = nG{x) [J^^G{x')dx']"~^ , where G is the unit 
Gaussian. For n > 1 both the mode and mean of P{x) are 
> 0. 

For an RTT the situation is more complicated, primar- 
ily because the numbers from which the extrema are chosen 
are correlated. In order to calculate an expectation value for 
the significance level of the extrema one would thus have 
to specify the covariance matrix which is exactly what we 
wanted to avoid. However, the expectation values of the ex- 
trema and their correlation are easily obtained from simu- 
lations. 

We thus performed 1000 simulations of the dataset {— 
10* spectra, hereafter SI) by randomly pla cing absorption 



o 
X o 



lines according to the parameters of Section 3.2 and using a 
constant S/N of 20. No proximity effect was included in the 
simulations. For each dataset we then constructed its RTT 
(using J — 00) and found the most significant deviations 
of the model from the data. The mean significance levels of 
these 1000 maxima and minima are shown as the dashed 
lines in Fig. ^ and the ±1(t regions are shown in grey. 

We determine the best fit value of J by calculating in 
the top panel of Fig. ^for each of the data points in the main 
panel. The best fit is achieved for log J21 = O.I^q *, where 
the errors are the formal 90 per cent confidence limits (outer 
dotted lines in Fig. ^). With Xmin ~ 1-98 the fit is acceptable 
{P{^ Xmin) = 0.16, where P is the x^-distribution with one 
degree of freedom). 

In order to check the error bars, the quality of the fit 
and the validity of the procedure as a whole, we performed 
a second set of 1000 simulations (S2). This time we included 
the proximity effect according to the ionization model of 
Section 3.3 with log J21 = 0.1. Each of these datasets was 
analysed in the same manner as the real data, i.e. for each 
we constructed Fig. ^ and measured log J. On average, the 
presence of the proximity effect is detected at the 5.2(tt 
level. The mean of the 1000 log J measurements is 0.10, 65 
per cent of the values lie within the range —0.15 < log J21 < 
0.35 and 86 per cent of the values lie within the range —0.3 < 
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Figure 3. Main panel: most significant positive (dots) and neg- 
ative (diamonds) deviations (in units of ax) of the absorption 
model from the data in the region most likely affected by the 
proximity of the QSOs as a function of the mean Lyman limit 
background intensity J. The dashed lines and the grey regions 
mark the expected and ±lcr significance levels of these deviations 
in the case where they are due to random fluctuations (i.e. for 
the correct model), which were determined from 1000 simulated 
datasets. The absence of a proximity effect corresponds to J = 00 
and is strongly rejected. Top panel: of the corresponding points 
in the main panel. The best fit is achieved for log J21 =0.1 with 
■ =1.98. The inner and outer dotted lines in the main panel 
mark the formal 68 and 90 per cent confidence intervals respec- 
tively. 



log J21 < 0.5. Finally, the fraction of measurements with 
Xmin ^ 1-98 is 0.10. We thus conclude that our method 
works well and that the error and quality of fit estimates 
above are reliable. 



4.2.1 Dependence on model parameters 

Does the above result depend sensitively on any of the model 



parameters introduced in Sections and 3.3? Because of 
the transition ( p^ we must expect the result to depend on 
/3. At a given J, a larger value of /3 increases the model 
transmission and thus decreases both the maxima and min- 
ima of Fig. ^ which will result in a larger measured value 
for J. For /? = 1.7 and 1.3 we find log J21 = 0.6+0:25 and 
— 0.451q 4^ respectively. However, increasing f3 also has the 
effect of decreasing the model variance so that the maxima 
and minima of Fig. ^ move further apart, which decreases 
the goodness of fit. For (3 — 1.7 we find x^ = 10.1 and thus 
the fit is no longer acceptable. 

Since the redshift coverage of the present sample is not 
very large and because we determine the optical depth nor- 
malization, B, directly from the data, our results cannot 
depend sensitively on 7. To confirm this we repeated the 
above analysis for 7 = 2.5 and found log J21 = 0.15lo3^. 

Cosmological parameters enter the analysis via the last 
factor in equation (|l2|). h cancels out but the denominator 
is relatively less sensitive to go than the numerator. For an 
open Universe this factor is larger than in the fiat case and so 
a larger J will be measured. For go = 0.15 we find log J21 = 
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Figure 4. Estimated value of the mean Lyman limit background 
intensity, J, assuming that the measured high ionization emission 
line redshift of every QSO is smaller than the true systemic value 
by an amount Av. Dots, solid line: this work. Diamonds, dashed 
line (offset by 50 km for clarity): Espey's (1993) analysis of 
Lu et al.'s (1991) data for comparison. Error bars are 90 per cent 
confidence limits. 



O.lSlEJ-g and thus our result does not depend sensitively on 



the cosmological model. 



4-2.2 Emission line shifts 

In Section ^ we already noted that the redshifts of the QSOs 
have probably been underestimated since they were deter- 
mined from high ionization lines. Assuming that all the emis- 
sion line redshifts are offset from their true systemic values 
by a velocity At) to the blue we have repeated the above 
analysis as a function of A«. In Fig. ^ we show our estimate 
of J for various values of Aw (dots, soli d line). The r educ- 
tion of J is very similar to that found by Espey (1993 1 (dia- 



mond s, dashed line). From similarly luminous Q SOs 



(1993) estimated t he mean velocity shift for the 



Espey 



Lu, Wolfe 



h Turnshek (199l| ) QSOs to lie in the range 1300 km s~^ 
< Au < 1600 km s~^ and estimated the true background 
intensity to be log J21 = — 0.31q'22- Fo'" ou r sample of QSOs 



the Ai;-luminos ity relationship given by Cooke, Espey, & 



Carswell (1997) predicts Aw ~ 800 km s~ which yields 
log J21 = —0.45^0 3- These results are in good agreement 
with each other as well as with log J21 ~ —0.3 a,t z = 2.5 



3 computed by Haardt & Madau (1996) for a background 



dominated by the QSO population observed in optical sur- 
veys and a go = 0.1 cosmology. 

4.3 Variation with luminosity and redshift 

4-3.1 Significance 

In Section we noted that the significance of the proxim- 
ity effect varied somewhat when excluding individual spec- 
tra from the combined RTT. We now examine the proximity 
effect in individual spectra in order to test whether it is cor- 
related with QSO Lyman limit luminosity or redshift. The 
former correlation would be expected if the proximity ef- 
fect were due to the increased ionizing flux in the vicinity 
of QSOs and the latter if in addition J varied with redshift 



and/or the luminosities and redshifts of the QSOs were cor- 
related. 

We thus constructed the RTTs of the individual spectra 
(using J = 00) and searched for the most significant posi- 
tive deviations of the m odel from the data in the same way 
as described in Section 4.2. In other words, for each QSO 
we determined the significance of the proximity effect. In 
Fig. |^(a) we plot these maxima as solid dots against QSO 



Lyman limit luminosity (calculated from equation ( |ll[ ) and 
listed in Table ^). For comparison, we analysed two sim- 
ulated datasets, SI and S3, in exactly the same way. SI 
was already introduced in Section (no proximity effect 
model). S3 incorporates the proximity effect according to the 
ionization model of Section 3.5 with log J21 = —0.45 and as- 
sumes that all QSO redshifts have been underestimated by 
Av = 800 km s~\ Writing M = Ma.x[RToo], we plot as open 
circles the mean of M found in the 1000 simulated datasets 
that include the proximity effect (S3). The open squares are 
the same for the no proximity effect simulations (SI). 

Previous authors have presented similar plots (e.g. 
EDO's Fig. 1) where they plotted the relative deficiency of 
absorption lines within some constant radius of the QSOs. 
One of the advantages of our method is that this radius is 
no longer constant but is rather allowed to vary in order to 
maximise the significance of the missing absorption. 

Let us first examine correlations in the simulated data. 
For SI there is clearly no correlation with luminosity which 
is as it should be. For S3 there seems to be a trend of in- 
creasing significance with increasing luminosity. However, 
the correlation does not seem to be as tight as one might 
expect and there seems to be no well-defined relation be- 
tween the two. Considering that these points are the mean 
of 1000 simulations this can hardly be due to random error. 
The only other parameter that varies from QSO to QSO in 
the simulations is the redshift. In Fig. ^(b) we plot the sig- 
nificance levels against QSO redshift. Again, for SI there is 
no trend. S3 exhibits the same sort of loose correlation as 
in Fig. ^a) . The luminosities and redshifts of the QSOs are 
not significantly correlated. Recall also that J does not vary 
with redshift in our simulations. Why then should there be 
a trend with redshift at all? The reason is that a given un- 
derdensity of absorption is more significantly detected when 
the 'background' absorption line density is higher than when 
it is lower (cf. LWC, Fig. 6b). Thus a given QSO will have 
a more noticeable proximity effect at high redshift (where 
the line density is higher) than at low redshift, all else being 
equal. 

We thus surmise that the lack of a well-defined relation 
between the significance of the proximity effect, M, and lu- 
minosity for S3 in Fig. ^(a) is due to the variation in redshift 
of the QSOs and that M is a function of both luminosity and 
redshift even though the QSO luminosities and redshifts are 
not significantly correlated and J does not vary with red- 
shift. 

To demonstrate this behaviour we first fit M with the 
function 



f{Lv,ZQ) = c^r itilL^ + m^ZQ 



(22) 



and then scale the significance levels to the mean redshift, 
z, and mean luminosity, L^. The results are plotted in 
Figs P(c) and (d) respectively. In Fig. |^ we plot M against 
ruLLv + mzZQ. We can now see that S3 exhibits almost per- 
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Figure 5. (a) Significance of proximity effect in individual QSOs versus QSO Lyman Umit luminosity. Solid dots are the data. Open 
circles are the mean significance levels found in 1000 simulated datasets, where the simulations include a proximity effect according to 
the ionization model with log J21 = —0.45 (S3). The vertical error bars on the data are the ±lcr ranges found in these simulations. Open 
squares: same as open circles but with J = 00 (no proximity effect, SI), (b) Same as (a) but now plotted against QSO redshift. (c) Same 
as (a) but all significance levels have been scaled to the mean redshift of the sample using equation (^) in order to isolate the effect 
of luminosity on the proximity effect. The solid, dashed and dotted lines show the best fit /(L^,5q) for the data and the two sets of 
simulations respectively, (d) Same as (b) but all significance levels have been scaled to the mean luminosity of the sample using equation 
( |22[ ) in order to isolate the effect of redshift on the proximity effect. The solid, dashed and dotted lines show the best fit f{Lv, zq) for 
the data and the two sets of simulations respectively. 



feet correlation with both luminosity and redshift and that 
the linear model gives a reasonably good description of 
the simulations. 

The discussion above implies that in order to properly 
disentangle possible correlations of the proximity effect with 
luminosity and redshift they should be determined jointly, 
not separately. We now investigate this effect in more detail. 

First, we need to choose a correlation statistic. Previous 
authors have often used Spearman's rank correlation coeffi- 
cient. However, in our case it is reasonable to assume that 
M is approximately Gaussian distributed and this is indeed 
observed in the simulations. Therefore it is not necessary to 
restrict ourselves to a non-parametric test. In addition we 
have already shown that the linear model (M) gives a good 
description of the simulations and from FigTq we can judge 
that it is also an acceptable model for the real data (see 
also Fig. ^). It is therefore reasonable to use the slopes of 



a linear fit as correlation measures. This also has the ad- 
vantage that the arguments of the following paragraphs can 
be understood analytically. In any case, we found that these 
arguments are qualitatively reproduced when using Spear- 
man's rank correlation coefficient. 

The simplest thing we can now do is to determine the 
slopes for the M-L^ and M-zq relations independently by 
fitting the functions 

to M. We have done this for the 1000 simulated datasets of 
both SI and S3 as well as for the real data. The results are 
shown in Fig. ^(a). The shaded regions in the main panel 
are the 68 and 95 per cent confidence regions for SI, the 
contours are the same for S3 and the cross marks the result 
for the real data. 
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Figure 7. (a) The dark and light grey shaded regions show the 68 and 95 per cent confidence regions of (m^ , 



) derived by independently 



fitting the significance of the proximity effect in the 1000 simulated datasets of SI (no proximity effect) with (t3)- The contour lines are 
the same for S3 (ionization model) . The cross marks the measured values and formal Iff errorbars of (m^ , m^) for the real data. The 
dashed lines correspond to the lines mriz = and = in panel (b). The upper and left panels show the one-dimensional probability 
distributions of m'^ and respectively. The dark and light shaded regions in these panels are the one-dimensional 68 and 95 per cent 
confidence regions respectively, (b) Same as (a) for (mz,mL), derived by fitting the significance of the proximity effect with 




"'L'-j/ ^ "'z^Q 

Figure 6. Significance of proximity effect in individual QSOs 
versus a linear combination of QSO Lyman limit lurninosity and 
redshift. Symbols have the same meaning as in Fig. H. The solid 
and dashed lines show the best fit /(L^,2;q) for the data and S3 
respectively. SI is not shown because for these points the best fit 
slopes are m/, niz ~ 0. 



The first thing we notice is that the values of (m^, m^) 
measured in the real data are consistent with the ionization 
model, but they are inconsistent with the no proximity ef- 
fect model at > 3a. However, note that in the simulations 

and m'^ are not independent. Calculating Pearson's cor- 
relation coefficient we find 



p(m'L,m'z) = piL^, zq) = 0.14, 



(24) 



where p{L^, zq) is the correlation coefficient between L^, and 
ZQ for the 10 QSOs used here, not that of the parent popula- 
tion (= p{Lij, Zq)) from which they were drawn, which may 
well be zero. Thus any measurement in the top right-hand 



part of the plot deviates from the no proximity effect hypoth- 
esis less significantly than what would have been inferred if 
the correlation between and zq had been neglected. Note 
that it is the actual numerical value of p that matters here 
and not whether or not it is consistent with p{L^, zq) = 0. 

What would we expect if the proximity effect were 
caused by some property, x, of the QSO or its environment 
unrelated to the QSO's UV flux ('a^-model')? As for the ion- 
ization model we would expect a correlation of AI with zq 
because this is simply the effect of increasing absorption 
line density with increasing redshift. Since p{L^, zq) 7^ we 
therefore expect the centre of the shaded ellipses in Fig. ^(a) 
to move along the line 



tu'l =mlm'z = p{L^, zq) 



V(L,) 



(25) 



(marked by a dashed lined), where V denotes the sample 
variance. Thus a measurement in the top right-hand part 
of the plot deviates less significantly from the a;-model than 
what would have been inferred if the correlation between Lv 
and ZQ had been neglected. 

This complication can be avoided if we use {niz , rriL ) 
of equation ( p^ instead of {m'z,m'j^). Fig. ^b) shows the 
result of fitting (^2|) to the simulations and the data. The 
relationship between Figs ^(a) and (b) can be most easily 
understood by writing down the linear co-ordinate transfor- 
mation which relates (m'j.,m^) and (mz,m_L) as 



ruz 
rriL 



l-p2 



-niz 
1 



(26) 



We can now see that they are related by a Lorentz transfor- 
mation followed by a stretch of (1 — p'^)" 2 and that the line 
m'j^ — m'z is transformed to the line m_t ~ 0. Thus in 
this frame of reference we do not have to worry about p when 
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asking whether the data is compatible with the x-model, ex- 
cept for the fact that we now have p{mL,mz) — —p{Lv, zq). 

In Figs^a), (b) and (c) we plot the distribution of min- 
imum values obtained from fitting ( p3[ ) and (|2^) respec- 
tively. The shaded histogram (SI) follows the x^-distribution 
quite well in all cases. However, both the real data (dotted 
line) and the ionization model (S3, solid histogram) are not 
well modelled by the independent fits (p3|). 

Thus we conclude that the observed significance of the 
proximity effect is a linear function of both redshift and 
luminosity and is well described by a function of the form 
(p^. This observed correlation with luminosity anrf redshift 
is inconsistent with the no proximity effect model at the 
3.5a level. A model that exhibits correlation with redshift 
but not luminosity is excluded at the 2.6<j level. These values 
take into account the correlation between and zq of the 
present dataset. If we had ignored this correlation we would 
have inferred 4.1(t and 3.2a respectively. 

The discussion above enables us to go one step further. 
Consider again the x-model. If x is indeed uncorrelated with 
for the general QSO population, i.e. p(x,Lv) = 0, then 
for our sample we would most likely find p{x, L^) = 0, which 
we have implicitly assumed in the previous discussion. How- 
ever, p{x,Lv) may well be 7^ 0, either because of random 
fluctuations or because p{x,L^) 7^ 0. This would induce 
a spurious correlation between M and L^. What value of 
p{x,L^) is required so that the x-model is consistent with 
the data? Assuming p(x,zq) = we find that p{x,Lv) has 
to lie in the range 



0.37 < p{x,L^) < 0.80 



(27) 



in order to be consistent with the data at the la level. Thus 
in the present dataset, the hypothetical property x would 
have to be noticeably correlated with luminosity. 

Since we do not know the distribution from which x is 
drawn it is not possible to reliably estimate whether this 
result excludes the hypothesis that x and L^, are not cor- 
related in the parent population. However, assuming Gaus- 
sianity in both x and L^, we find that the most likely value 
of p{x,L^) — 0.59 excludes the hypothesis p{x,L^) — 
at the ~ 96 per cent confidence level (using Student's t- 
distribution) . 

By comparing the data with S3 in Fig. 0(b) it is ap- 
parent that the observed dependence of the proximity effect 
on redshift is entirely accounted for by the evolution of the 
number density of absorption lines and there is no evidence 
to suggest that J varies over the redshift range considered 
here. It is also apparent that this conclusion is not overly 
sensitive to our choice of 7, the evolutionary index of the 
absorption line density. On the other hand, the data do not 
exclude some (positive or negative) evolution either. In any 
case, if J is due to the known QSO population then it is ex- 
pected to peak smoothly in the redshift range covered here 
( Haardt fc Madau 199(\ ) and thus one would not expect 
strong evolution. 



4.3.2 Size 

If the ionization model is correct one would in principle also 
expect the size of the region affected by the QSO's UV flux to 
be correlated with luminosity. We can test for this effect by 




Figure 8. (a) The grey shaded histogram is the distribution of 
minimum values obtained from fitting f]^{Lu) (equation pj]) 
to the 1000 simulated datasets of SI (no proximity effect). The 
solid histogram is the same for S3 (ionization model). The dashed 
line marks the value of x^ for the real data. The smooth solid line 
shows the x^-distribution for eight degrees of freedom, (b) Same 
as (a) for fz{zQ). (c) Same as (a) and (b) for f{L^, zq) (equation 
^2! ). The solid line shows the x'^-distribution for seven degrees of 
freedom. 



noting the size of the smoothing function at which the max- 
ima M of Fig. P were detected. However, in the simulations 
we flnd that the distribution of sizes is very broad and almost 
one-tailed and thus difficult to characterise. In addition the 
distributions for SI and S3 overlap almost completely with 
only the peak moving to slightly larger values for S3. There 
is also some evidence that for the ionization model the distri- 
bution of sizes moves to larger values for larger luminosities 
but again the effect is small compared to the extent of the 
distribution. Thus we are forced to conclude that this is not 
a powerful test and simply note that the detected sizes in 
the real data range from ~ 500 to ~ 7000 km s~^ and lie 
well within the 90 per cent confidence ranges of both SI and 
S3 in all cases. 



5 RESULTS ON FOREGROUND PROXIMITY 
EFFECT 

5.1 Existence of the effect 

We now exploit the fact that the QSOs of Table |l| are a 
close group in the plane of the sky. Essentially, we repeat 
here the analysis of Section 4.1: for each pair of background 



QSO (BQSO) and foreground ionizing source (IS), we shift 
the spectrum of the BQSO into the restframe of the IS 
and construct its transmission triangle using J = 00. We 
now use a Gaussian smoothing function because we are no 
longer working at the 'edge' of the data. We then average 
all the transmission triangles where they overlap, creating a 
composite RTT. Thus in equation (EG) Ti still refers to the 
measured transmission of the BQSO(but is replaced by d 
because we now use a Gaussian smoothing function) and we 
replace zqi with zisi, where i now labels a BQSO-IS pair. 
If the foreground proximity effect exists, it should be more 
significant in this composite RTT than in any individual 
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log cj^ 

Figure 9. Significance of tlie foreground proximity efi'ect as de- 
rived from composite RTTs which include all BQSO-IS pairs with 
u)p > oJc- The solid line is the data. The dotted line and grey 
shaded regions are the mean and ±lcr significance levels found in 
1000 simulated datasets that include the foreground proximity ef- 
fect according to the ionization model with log J21 = —0.45 (S4). 
The dashed line and hashed region are the same for J = 00 (no 
proximity effect, SI). The number of BQSO-IS pairs included in 
the RTT increases from 1 at the left end of the plot to 14 at the 
right end. 



RTTs and is expected to appear as a region of RGoc > 
near Ar = A^. From now on we exclude all spectral regions 
within 5000 km s~^ of the BQSOs from the analysis in order 
not to contaminate our results with the background prox- 
imity effect. 

Before we proceed we need to consider the follow- 
ing complication. For each BQSO-IS pair we can calculate 
LUp = uj{z = 2:is,2:is,a) (assuming log J21 = —0.45, cf. equa- 
tion ^2]). ijp is the maximum value a given IS can achieve 
along the line of sight to a BQSO, separated on the sky by 
an angle a. Pairs with large cup should show a strong prox- 
imity effect. However, below some value of uip the proximity 
effect will be essentially non-existent. Adding pairs with Ljp 
below this value to the composite may in fact decrease the 
overall significance of the effect. 

We thus construct several composite RTTs, each with 
a different lower limit on cjp, which we denote by luc- Thus 
RGj{luc) refers to a composite reduced transmission triangle 
in which all BQSO-IS pairs with ujp > Uc are included. 

Each of these RTTs was searched for the most signif- 
icant positive deviation of the model from the data in the 
interval [-1000 km s"\ +2000 km s"^] around A^,. The in- 
terval is asymmetric to allow for underestimated QSO red- 
shifts. In Fig. ^ we plot these maxima versus uic- For com- 
parison we have performed the same analysis on the 1000 
simulated datasets of SI (no proximity effect). We plot the 
mean and dzla significance levels as the dashed line and 
hashed region respectively. Clearly, our data are consistent 
with the absence of any foreground proximity effect. 

Given the luminosities and inter-sightline spacings of 
the present set of QSOs, do we expect to be able to detect 
a signal? To answer this question we have created a fourth 
set of 1000 simulations (S4). These include the effects of all 
foreground IS with ujp > 0.5 and log J21 = —0.45. Subjecting 
S4 to the same analysis as above yields the dotted line and 
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Figure 10. This plot is the equivalent of Fig. ^ for the foreground 
proximity effect. In the main panel we show the most significant 
positive (dots) and negative (diamonds) deviations (in units of 
ctq) of the absorption model from the data in the regions near 
the positions of close-by foreground QSOs as a function of the 
mean Lyman limit background intensity J. The dashed lines and 
the grey regions mark the expected and ztlcr significance levels of 
these deviations in the case where they are due to random fiuctu- 
ations. The data are consistent with the absence of a foreground 
proximity effect (J = 00). Thus we cannot derive an upper limit 
on J. The two dotted lines mark the 68 and 90 per cent confidenc 
lower limits on J. 



grey shaded region in Fig. ^. Evidently, if J has indeed the 
value that was measured from the background proximity 
effect and if QSOs radiate isotropically, then we should be 
able to detect the foreground proximity effect at the 2 — 4cr 
level in our data. 



5.2 Lower limit on J 

The above result implies that we can at least derive a lower 
limit to J under the assumption that QSOs radiate isotropi- 
cally . Se tting ujc = 0.5 we have repeated the analysis of Sec- 
using a foreground RTT that includes 14 BQSO-IS 



tion 4.2 ^ ^ 

pairs. The result is shown in Fig. |l^. We can see that the 
data are consistent with a large range of values of J, includ- 
ing J = 00. However, for small values of J the model predicts 
too little absorption to be compatible with the data. From 
this constraint we derive a lower limit of log J21 > 0.3 (90 
per cent confidence). 

From Fig. ^ we can see that this lower limit is larger 
than the upper limit derived from the background proxim- 
ity effect for all At; ^ 200 km s'V If we did not underes- 



timate J in Section 4.2 then the simplest explanation for 
this discrepancy is that QSOs radiate anisotropically. If we 
require that the lower limit derived from the lack of a fore- 
ground proximity effect coincides with the upper limit of 
the background measurement then the QSOs of our sam- 
ple must emit less ionizing radiation in the plane of the sky 
than along the line of sight to Earth by at least a factor of 
1.4 for Av > 400 km s~^. This number increases to 2.2 for 
Av > 800 km s"\ 
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Table 2. The Q0042-2639 quadrangle. 



BQSO 



a M 
arcmin ctq 



A FWHMs 
km s"-*^ km 



Q0043-2633 


13.6 


2.97 


-460 


1260 


Q0042-2627 


11.7 


3.11 


-2390 


5530 


Q0041-2638 


11.9 


2.71 


+230 


3120 


Q0042-2656 


17.6 


3.03 


+2630 


1210 


Composite 




3.52 


+160 


1880 



5.3 The Q0042-2639 quadrangle 

Having established that the dataset as a whole does not ex- 
hibit any evidence for the existence of the foreground prox- 
imity effect, we now present a possible exception to this rule. 
From Fig. □ we can see that Q0042-2639 is surrounded by 
four nearby background QSOs. Interestingly, all four show 
underdense absorption near the position of the foreground 
QSO. In Table ^ we list the BQSO name, the angular sep- 
aration from the foreground QSO, the significance of the 
underdensity, M (in units of (Jg), the velocity offset of the 
underdensity from the foreground QSO redshift and finally 
the size of the underdensity (FWHM of the smoothing Gaus- 
sian). In the last line we list the same quantities for the 
composite RTT of the four BQSOs. Note that the offsets 
from the foreground QSO's redshift are of the same mag- 
nitude but of opposite sign for the two pairs of opposing 
BQSOs (cf. Fig. p|). However, this is no longer true if we add 
^ 400 km s^^ to the foreground QSO's redshift. 

At first glance, it may seem exceedingly unlikely to find 
an underdensity in four different lines of sight at a simi- 
lar position (which happens to coincide with a foreground 
QSO) by chance. However, the spread of the underdensities 
in redshift is actually fairly large. Thus when we combine 
the four lines of sight to a composite RTT the significance 
level of the 'void' rises only marginally from ~ Sao in indi- 
vidual lines of sight to 3.5aG in the composite. Nevertheless, 
in the 1000 simulated datasets of SI we find that the largest 
random fluctuations in the composite RTT have a mean sig- 
nificance level of (1.69 ± 0.53)(tg and thus it seems that the 
void is in fact genuine. Indeed, a measurement of J from the 
composite RTT yields log J21 = — CGSlg g which is lower 
than, but consistent with our measurements from the back- 
ground proximity effect. From Fig. ^ (and Table |l]) we can 
see that Q0042-2639 also shows a slightly (but not signif- 
icantly) more prominent background proximity effect than 
'predicted' by the other QSOs for its luminosity and redshift. 



Together these observations may be an indication that the 
luminosity ot Q0042-2t)il9 has been underestimated. 



6 UNCERTAINTIES 

In Section 



4.2 we have already quantified the effect of uncer- 
tainties in model parameters, QSO redshifts and cosmology. 
We now discuss a number of other uncertainties associated 
with the present study. 

It is clear that our method is sensitive to errors in the 
continuum placement, pro bably more so th an the generic 
'line counting' method. In Liske et al. (200C| ) we attempted 
to account for random errors in the lowest orders of the con- 



tinuum fit by determining the normalization of the mean op- 
tical depth, B, for each spectrum separately. Here we used 
a single value of B derived from the entire dataset. This has 
the advantage of making our estimate of J fairly insensitive 
to the adopted value of 7 (the evolutionary power law index) . 
The disadvantage is that errors in the continuum placement 
should, in principle, increase the scatter when comparing the 
significance of the proximity effect in individual lines of sight 
(cf. Fig. However, this comparison is based on measure- 
ments near the Lya emission lines of the QSOs where the 
S/N is in general quite high and thus the continuum more 
secure. A systematic over- or underestimation of continua 
can only affect our results if such a bias is different for dif- 
ferent parts of the spectra. This difference could be caused 
by the higher S/N and greater curvature of the continuum 
in the wing of the Lya emission line. 

As our method deals only with the transmitted flux it 
circumvents all problems that arise from deflning a sample 
of individual absorption lines. These problems include line 
blending, curve o f growth effects (e.g. ^c ott_et_al . 2000 ) and 
Malmquist bias (Cooke, Espey, & Carswell 1997). 



The environment of QSOs may well differ from the in- 
tergalactic environment in other aspects than just the in- 
tensity of ionizing radiation. Most importantly, there may 
be additional absorption in the vicinity of QSOs above and 
beyond the absorption already accounted for in our model. 
If QSOs are hosted by groups or clusters of galaxies then 
the gravitational pull of the host will cause infall of the sur- 
rounding material an d may thus increase the a bsorption line 
density near a QSO (Loeb & Eisenstein 1995). If this effect 



is not taken into account, the proximity effect will appear 
weaker than it should and a larger value of J will be inferred. 
Loeb fc Eisenstein (1995[ ) suggested that the magnitude of 
this effect may be as large as a factor of ~ 3. If QSO luminos- 
ity is correlated with host mass, then more luminous QSOs 
are affected more strongly by clustering and the value of J 
derived from the brightest QSOs should be higher than that 
derived from the faintest. On the other hand, if all QSOs 
were affected by clustering in the same way, then one might 
expect the observed slope of the significance-luminosity rela- 
tion in Fig. ^ to be larger than the one expected on the basis 
of the measured (and overestimated) J. This is actually the 
case, although not significantly so. However, all of the QSOs 
in the present sample are radio-quiet. The hypothesis that 
such QSOs reside in rich galaxy cluster environments has 
been repeatedly rejected at low redshifts and there seems to 
be little evolutio n in the environment of radio-quiet QSOs 
up to z < 1. 5 dCroom fc Shanks 1999|; |Smith, Boyle, & 



Maddox 2000). In any case, it is most likely very difficult 
to disentangle clustering from the proximity effect and we 
beheve that this issue deserves further study. For now, it 
remains an uncertainty. 

In Section 3.3 we assumed that the column density of 
an absorber is proportional to the inverse of the incident 
ionizing flux. This is onl y true for absorbers composed purely 
of hydrogen. However, scott et al. (2000) found that the 
inclusion of metals into the model has an insignificant effect 
on the derived value of J. 

We have also ignored the fact that as radiation travels 
from a QSO to a given absorber it will be attenuated by 
all the intervening absorbers. In particular, an intervening 
strong Lyman limit system will essentially 'black out' the 
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QSO entirely. In princip le, disregarding this effect cause s 



overestimation of J but Cooke, Espey, & Carswell (199'i) 
concluded that the effect is negligible for log TV < 17. Higher 
column densities produce a conspicuous continuum break at 
the Lyman limit. Unfortunately, only three of our spectra 
cover any part of the Lyman limit region. One of these, 
Q0042-2656, shows a Lyman break but the system lies ~ 
5200 km s~^ away from the QSO. In any case, since high 
column density absorbers are comparatively rare we would 
expect that only a fraction of our QSOs are afflicted by this 
problem. However, from Fig. |^ it is apparent that none of 
our QSOs shows a proximity effect which is unusually small 
for its luminosity and redshift. We thus find it unlikely that 
we have significantly overestimated J due to this effect. 

Obscuration by dust in intervening damped Lya ab- 
sorption systems may cause a QS O's luminosity and conse - 
quently J to be underestimated (srianand & Khare 1996). 
Since there are no known damped Lya absorption systems 
or candidates in our sample we believe that our value of J 
is not affected by dust obscuration. 

However, there are a number of other factors that cre- 
ate uncertainty in the estimated luminosities of the QSOs. 
Errors in 7^- corrections and the QSO continuum slope could 
be avoided by direct spectrophotometric observations, but 
there are at least two other more fundamental and probably 
larger uncertainties: 

1. QSO variability. The equilibration time-scale of Lya 
absorbers, icq, is of the order 10* years. Thus the observed 
ionization state of an absorber in the vicinity of a QSO will 
approximately reflect the ionizing flux received from that 
QSO averaged over the 10* years prior to the epoch of ob- 
servation, L. Let us assume that the intrinsic luminosity 
of QSOs varies on a single time-scale, tv. For tv ^ toq the 
QSO's observed luminosity may be different from L and thus 
the strength of the observed proximity effect may be differ- 
ent from that expected on the basis of the ionization model. 
The same may be true for the foreground proximity effect 
even when Iy > icq if iv is smaller than the light travel time 
from the IS to the observed absorption (« 10^ years). Ob- 
viously, we have no informa tion on QSO variab ility on such 
large time-scales. However, Giveon et al. (199£| ) found that 
on tim e-scales of 100 to 1000 days the distribution of bright- 



ness de dations about the mean light curves of 42 QSOs has a 
width of 0.14 mag in the B-band. Thus, variability on short 
time-scales contributes substantially to the uncertainty in 
the QSOs' luminosities. 

2. Gravitational lensing. None of the QSOs in our sam- 
ple are known to be lensed and at least four of them have 
been included in searches for mult iply imaged QSOs with 
negative results ( ^urdej et al. 199S ) , implying that they are 
at least not strongly lensed. However, some or all of them 
may be weakly lensed by the non-homogeneous distribution 
of foreground matter on large scales. There are two known 
overdensities of C iv absor ption systems in the foreground of 
the present QSO sample (Williger et al. 1996) which could 
possibly ca use weak lensing ([Vanden Berk, Quashnock , & 



j9|; [Holz fc Wald 19981 ). Ray-tracing experiments 



York 1 

performed on simulations of cosmological structure forma- 
tion have yielded the probability distribution function of 
the magni flcation caused by galaxy clusters and large-scale 



tribution can be as high as 0.4 (with a mean of 1) and the 
distribution has a power law tail towards large magnifica- 
tion. For the standard model the probability of encountering 
a magnification of 2 or greater is a few per cent but for a 
A model it is less than 10~^. A magnification by a factor of 
1.2 (or a demagniflcation by a factor 1.2~^) introduces an 
additional uncertainty of ±2.5 log 1.2 — ±0.2 mag which is 
actually larger than the quoted measurement errors on the 
observed B-band magnitudes. 

What are the effects of luminosity uncertainties on our 
results if we treat them as random errors? Any statistical 
errors should increase the error on J. We can estimate the 
magnitude of this effect by the following argument. If we 
measured J from a single QSO then the additional error 
on J should be on the order of A(log J) « 0.4Am, where 
Am is the error on the QSO's magnitude. For ten QSOs 
this error should be smaller by a factor of vTO. Thus for 
Am = 0.5 mag we find A (log J) — 0.063 which is much 
smaller than the quoted error of 0.4. 

This also shows that the known QSO variability, gravi- 
tational lensing and measurement errors on the QSO magni- 
tudes cannot increase the errorbars on our foreground and 
background estimates of J by an amount large enough to 
make these two values compatible. 

However, statistical errors in L,^ will weaken our results 
on the correlation between the significance of the proximity 
effect and L^. Clearly, for errors as large as Am = 0.5 mag 
(which corresponds to a factor of 1.6 in luminosity) fitting 
a straight line to the data points in Fig. |^ will be almost 
meaningless. In Section we estimated that the signifi- 
cance of the proximity effect is correlated with luminosity 
at the 2.6a level. For Am — 0.3 mag this significance drops 
to 2.1(7. Thus measurement errors and the known variabil- 
ity of QSOs alone cannot entirely invalidate the evidence for 
a correlation of the background proximity effect with QSO 
Lyman limit luminosity. 

In any case it is probably more appropriate to treat 
both QSO variability and gravitational lensing as system- 
atic errors. Because of the steepness of the bright end of 
the QSO luminosity function a magnitude limited sample 
of QSOs is more lik ely to contain magnified QSOs than de- 
magn ified ones (e.g. |^ei 1995 ; Hamana, Martel, & Futamase 



200C). The same is t rue for QSOs that are near a peak in 



structu re (Hamana, Martel, fc Futamase 2000 ; Wambsganss 



Cen, & Ostriker 1998|) . At z = 3 the dispersion of this dis 



their lightcurves (e.g. Francis 1996 ). Thus either of these ef- 
fects could cause a systematic overestimation of J. Together 
they imply that on average the QSO magnitudes may have 
been overestimated by ~ 0.35 mag which corresponds to an 
overestimation of J by A (log J) = 0.14. 

Note that gravitational lensing and QSO variability on 
short time-scales (< 10* years) affect our J estimates from 
the classical and foreground proximity effects in the same 
way. However, if QSO luminosities vary on time-scales of 
~ 10® years then this will affect only the latter estimate as 
explained above. In this case we would expect the QSOs to 
have been systematically fainter over a period of ~ 10® years 
prior to the time they emitted the photons which we receive 
today, causing us to overestimate J when measuring it from 
the foreground proximity effect. To reconcile the foreground 
and background values of J we require a variability of Am = 
2.5 log 2.2 — 0.86 mag (assuming At; = 800 km s~^) on time- 
scales of ~ 10® yearS; 



Past authors ( Cooke, Espey, fc Carswell 1997 ; Scott 
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et al. 2000) have attempted to test for the presence of gravi- 
tational lensing in their data: since high luminosity QSOs are 
more likely to be lensed than low luminosity ones, an esti- 
mate of J from the former group should be higher than from 
the latter. From Fig. |5| we can see that the higher luminosity 
QSOs of our sample actually show a slightly more prominent 
proximity effect than expected for log J21 = —0.45. Thus a J 
measurement from these four QSOs will yield a lower value, 
contrary to the expectation if they were lensed. 



7 CONCLUSIONS 

We have analysed the Lya forest spectra of a close group of 
10 QSOs in search of the (foreground) proximity effect using 
a novel method based on the statistics of the transmitted 
flux. We list our various measurements of J in Table ^ and 
we summarise our main results as follows: 

1. We confirm the existence of the classical background 
proximity effect at the > 99 per cent confidence level. 

2. From the observed underdensity of absorption near 
the background QSOs we derive log J21 = O.llg's (90 per 
cent confidence limits). 

3. Correcting all QSO redshifts by Aw = -1-800 km s"^ 
we find log J21 = — 0.45!to'3. The reduction of J with Av is 
consistent with previous results. 

4. The significance of the background proximity effect in 
individual lines of sight is correlated with QSO Lyman limit 
luminosity at the 2.6a level, thus lending further support 
to the hypothesis that the proximity effect is caused by the 
additional UV ffux from background QSOs. We account for 
the fact that the significance is also correlated with redshift 
which is due to the evolution of the absorption line density. 
We considered an alternative model for the proximity effect 
where the underdense absorption is caused by some hypo- 
thetical property x of the QSO or its environment. We find 
that the property x would have to be noticeably correlated 
with luminosity in our sample and probably in the QSO 
population in general. 

5. The full sample shows no evidence for the existence 
of the foreground proximity effect. 

6. This absence implies a lower limit of log J21 > 0.3. If 
we interpret the discrepancy of this lower limit with previ- 
ous measurements as evidence that QSOs radiate anisotropi- 
cally, then they must emit at least a factor of 1.4 less ionizing 
radiation in the plane of the sky than along the line of sight 
to Earth. 

7. Our sample includes the fortunate constellation of 
a foreground QSO surrounded on all sides by four back- 
ground QSOs with approximately equal separations from 
the foreground QSO. Contrary to the rest of the sample, 
this particular QSO induces a foreground proximity effect 
in the surrounding lines of sight at the 3.5a level. For this 
subsample we measure log J21 = — O.GSIq'j. 

8. Finally, we have discussed possible sources of sys- 
tematic and additional statistical errors. We conclude that 
clustering of absorption systems around QSOs is the most 
likely source of systematic error. We also find that the known 
variability of QSOs reduces the significance of the proximity 
effect-luminosity correlation to 2.1cr. 

From Fig. |l| we can see that our measurement of J 
from the background proximity effect is consistent with 



Table 3. Summary of J measurements and limits. 





log J21 


Errors" 


standard 


0.10 


-1-0.40, -0.30 


/3 = 1.3 


-0.45 


-1-0.45, -0.40 


13 = 1.7 


0.60 


-1-0.30, -0.25 


7 = 2.5 


0.15 


-1-0.35, -0.30 


go = 0.15 


0.15 


-1-0.40, -0.30 


Av = 400 km s'^ 


-0.25 


-1-0.40, -0.30 


Av = 800 km s'^ 


-0.45 


-1-0.40, -0.30 


foreground 


> 0.30 




Q0042-2639 quadrangle 


-0.65 


-1-0.60, -0.50 



"Errors are 90 per cent confidence limits. 



most previous measurements. Bechtold (1994 ) found a some- 
what larger value log J21 = 0.5 but she did not correct 
the QSO redshifts, remarking only that a correction of 
Av = 1 000 km s~^ would lower her value of J by a fac- 
tor of 3. Srianand & Khare (1996) also found a considerably 



higher va lue which is part ly due t o the fact that they used 
f3 = 1.7. [Kulkarni fc Faii| 's ( |l993[ ) low value of J21 = -2.2 
at z ~ 0.5 is usually taken as evidence for an evolving back- 
ground. 



Bechtold (1994) and Scott et al. (200C ) noted that the 



relative deficit of absorption lines near background QSOs 
was larger for the high-luminosity halves of their samples 
than for the low- luminosity ones. However, they did not 
quantify this effect in any detail and did not compare it to 
the expectations from the ionization model. BDO identified 
a trend of the line deficit with luminosity at la significance. 
Lu, Wolfe, fc Turnshck_(199lJ)_ on the other hand found no 
such trend. In Figs. ^, |^ and ^ we have presented good evi- 
dence that the significance of the proximity effect is indeed 
correlated with QSO luminosity and redshift. This result 
was possible because we used a technique that is more sen- 
sitive to variations of the absorption density on large scales 
than the generic line counting method. Using simulations 
we find that these correlations are entirely consistent with 
the expectations of the ionization model. However, consid- 
ering the discussion of Section |^ we clearly need to apply 
our method to a larger sample with a wider range in lu- 
minosity to establish this result more firmly. Nevertheless, 
the present analysis provides further evidence that the inter- 
pretation of the proximity effect as being due to increased 
ionization caused by the extra UV flux from the background 
QSO is essentially correct. 

Although their results were poorly constrai ned, it is in- 



teresting to note that the three BQSO-IS pairs of Fernandez- 



Soto et al. (1995) favoured a similarly high value of J as our 



full sample of 14 pai rs. The non-d etecti ons of the foreground 



proximity effect by Crotts (1989) and M0ller fc Kjffirgaard 
(1992| ) also imply a high value of J. On the other hand 



all possibl e positive detections of the foreground prq xim- 
ity effect ( [Pobrzycki fc Bechtold 199lt |Srianand 1997|; the 



Q0042-2639 quadrangle) yield J values that are in line with 
the measurements from the background effect. 

One possible explanation is that most of the QSOs were 
substantially fainter over a period of ~ 10® years prior to the 
time they emitted the photons which we receive today. This 
would cause an overestimate of J when measured from the 
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foreground proximity effect but would not affect the results 
of the background proximity effect. 

As we have already pointed out, the other explanation 
is that QSOs radiate anisotropically. There is a large body 
of observational evidence which suggests that the observed 
characteristics of an Active Galactic Nucleus (AG N) depend 



on t he d irection from which it is viewed (see e.g. Antonucci 
1993 | fo| r a review). The basic theme of unified models for 
AGN is that some or even all of the many different types of 
AGN are in fact the same type of object but seen from dif- 
ferent directions. In these models the directionality is caused 
by a thick, dusty and opaque torus which surrounds the cen- 
tral engine. The smooth continuum and broad emission lines 
of a QSO are thought to originate from within the torus and 
thus they can only be seen indirectly by scattered light when 
the torus is viewed approximately edge on. 

This scenario has two obvious implications: i) all QSOs 
show a background proximity effect and ii) whether a fore- 
ground proximity effect is seen or not depends on whether 
nearby absorption systems probed by other sightlines can 
'see' inside the torus of the QSO. If so, they will roughly 
see the same continuum as we do, resulting in a measurable 
depletion of absorption and a corresponding J value which 
is similar to that derived from the background effect. If not, 
there will be little or no foreground effect, resulting in a high 

value or lower limit for J. 

Assuming a simple picture of this kind, Dobrzycki & 



Bechtold (1991) used the velocity offset of their void from 
the foreground QSO to derive a value of ~ 140° for the 
opening angle of the torus. 

However, it is difficult to explain the case of the Q0042- 
2639 quadrangle with this scenario. If all of the four under- 
densities are real and caused by the forground QSO then 
it has to emit similar amounts of radiation along three 
nearly perpendicular axes (north-south, east-west and to- 
wards Earth), leaving little room for anisotropic emission. 

The current situation is thus uncertain and intriguing 
enough to stimulate further observations. The motivations 
and potential gain are clear: by investigating the radiative 
effects of QSOs (or other AGN) on nearby absorption sys- 
tems along other lines of sight we can 'view' them from dif- 
ferent directions which may help to constrain unified models 
of AGN. The ideal targets for further observations are dense 
groups of QSOs at similar redshifts. These are rare in current 
QSO catalogues but, fortunately, currently ongoing redshift 
surveys like 2QZ will soon remedy this situation. 
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